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Abstract 

The reactive process of barrier escaping from the metastable potential well is studied together 
with the extension of Kramers’ rate formula to the fractional case. Characteristic quantities are 
computed for an thimbleful of insight into the near barrier escaping and recrossing dynamics. 
Where the stationary transmission coefficient is revealed to be larger than the usual cases which 
implies less barrier recrossing. And the non-monotonic varying of it reveals a close dependence to 
the fractional exponent a. In most cases, the near barrier behavior of the escaping dynamics is 
equivalent to the diffusion in the two-dimensional non-Ohmic damping system. 
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I. INTRODUCTION 


Among all kinds of stochastic processes that produce sub-diffusion, fractional Brownian 


motion (fBm) may be the model which is most particularly relevant to the study of 
disciplinary sciences. For example, in the biochemistry research of subcellular transport, 
negative and long-eange cotnelat.on hae been obeeeved m auMiffneing .nRNA naoleculea g 
RNA proteins, and chromosomal loci within E. coli cells [3[ when the fractional exponent 


H is below 1/2. A 


translocations jd, 5|, the dispersion of apoferritin proteins in crowded dextran solutions 


eanwhile, fBms can also be used in a similar way to describe unbiased 


and lipid molecules in lipid bilayers j^. For these reasons, study on the model of fBm and 




its variances has attracted considerable attention in recent years 


12 |. 


However, we noticed that, in the currently existing studies of fBms, particles are always 
conhned in a harmonically varying bounding state. The barrier escaping of fBm particles 
from a quasi-bounded metastable potential well is seldom considered. While from another 
point of view, barrier escape problems are of fundamental interest in physics and chemistry. 
A great amount of chemical events, such as chemical reactions, molecular diffusion, or col¬ 
lision of molecular systems, can be modeled W a single barrier escape process within the 


framework of standard Brownian motion 


15l |. Great progress has been witnessed in the 


study of such kind of problems since the foundation of Kramers remarkable reaction rate 
formula [^. But no extension has been made to the fBm case. Therefore in this paper, we 
report a recent study of us concerning on the fractional Kramers problems which may kill 
these two birds with one stone. 

The paper is organized as following: in SecHIl a briefly mathematical derivation concern- 
ing on the extension of the rate formula is presented by using of the reactive flux method 
jlTI. Il8l|. In Sec lIIIl the time-dependent transmission coefficient is computed by Laplacian 
solving the generalized Langevin equation of the system. The conclusions of this study are 
summarized in Sec IIVI where some prospect are also made for further discussions which may 
be concerned in the near future. 
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II. REACTIVE FLUX FORMALISM FOR THE RATE 


For the expression of the escape rate, the traditional starting point throngh the method 


of reactive flux is alternatively to do the following integration 


17 


20 | 


k{t) = 


m 




^+oo 


Qh y_c 


dx[\ 


dvoW {xo,vo)vox{xo = Xb,VQ;t), 


( 1 ) 


where Q = Wst{x)dx is the partition function of the system for reactions integrating 
over the distribution of the ground state Wst{x). h denotes the phase cell and m the mass of 
the particle. W{xq,vq) = exp[—F7(a;o, no)/A:BT] is assumed to be the equilibrium distribution 
and H{xo,vo) = \m{ujlxQ + Vq) is the Hamiltonian of the system. 

The function is a characteristic quantity computed corresponding to an en¬ 

semble of diffusing particles which are assumed to start from xq = Xb- Mathematically it 
can be resulted from the following integration on the reduced distribution function 


x{xo,Vo;t) = 


W (x, xo, Vo] t)dx = -erfc 


where 


lU(x,Xo,no;^) = 




exp 



is reduced from the joint probability density function (PDF) of the system 21| 


W{x,V,XQ,Vo]t) = 




( 2 ) 


( 3 ) 


( 4 ) 


27r|A(t)|V2 

with y{t) the vector [x — {x{t)),v — (n(t))] and A{t) the matrix of second moments in which 
= ([x — {x(t))Y) is the variance of x(t). 

In the spirit of reactive flux calculation, the initial conditions are assumed to be at the 
top of the barrier. The rate is accounted from an ensemble of reactive trajectories which 


start with identical initial conditions but experience different stochastic histories. T 
the dynamics are contained in the generalized fractional Langevin equation (GFLE) 
of the system 

mx + [ r]{t — t')x(t')dt' + dxU{x) = ^{t) 


len all 


22 


24| 


( 5 ) 


where ri{t) = rjat “/F(l — a) is the frictional kernel with fractional exponent 0 < a < 1 and 
rja the strength constant, ^{t) is the fractional Gaussian noise with {^(t)^(t')) = ksTrialt — 


t'\ “ satisfying the Kubo’s fluctuation-dissipation theorem 


25]. 
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Eq.Q is one of the alternative form of GFLE for fBm systems which can also be derived 


mathematically from the system-plus-reservoir model of harmonic oscillators 


26 


271. Where 


U (a:) is the external potential which is always assumed to be a metastable form in the study of 
barrier escaping dynamics. But the approximation of an inverse harmonic potential is always 
made in the neighbourhood of the saddle point. Supposing the temperature of the system 
is much less than the barrier height of the potential, a Boltzmann form of the stationary 
probability distribution Wst{x) will be satished in the neighbourhood of the potential well. 
Then all the quantities can be obtained from Laplacian solving the GFLE. 

In order to isolate th e dy namical corrections to the transition state theory (TST) rate 

kTST ^ k^^-V^lksT 
Qh 
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30| . it is convenient to dehne the transmission coefficient 


K{t) = 


m 


ksT 


^H-oo 




( 6 ) 


Substituting Eq.( 


into the above integration we obtain immediately 


^ ^ keTH^it)) 


(7) 


where H{t) = £“^[( 5 ^- 1 - 577 ( 5 ) — 0 ;^)“^] is namely the response function obtained from solving 
the GFLE and 77 ( 5 ) = is the Laplacian transformation of the friction kernel 77 (f). This 

expression for the fractional reactive index leads immediately to Kramers formula for the 
rate constant [l^ which physically describes the probability of a particle already escaped 
from the well to recross the barrier. 


III. FRACTIONAL TRANSMISSION COEFFICIENT 


In the numerical computations here and following, all the quantities are rescaled to di¬ 
mensionless ones for simplicity. Firstly in FigJTl we plot the time-dependent transmission 
coefficient K{t) at various fractional exponent a. From which we can see that, the stationary 
value of K,{t) (here it is dehned as Kst = lirrit^ooY)) is comparable in most cases to the 


results of two-dimensional non-Ohmic damping systems 20||. However it is always larger 


than the usual cases such as that of one-dimensional non-Ohmic damping . This reveals 
that the fBm particles already escaped from the potential well will have less probability to 
recross the barrier. 
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FIG. 1: Time dependent varying of n{t) at various a. Parameters in use are rja = 2.0, w = ksT = 
1.0, xq = —1.0 and vq = 2.0. 




FIG. 2: Time dependent varying of {x(t)) and at various a. 

In order to find the inherent reason for the occur of this physical scenario, we give an 
illustration on the time dependence of {x{t)) and cr^(t) in Figj2l Seen from it that the mean 
displacement {x(t)) diverges negatively at a. However, the value of Anit) is always positive 
increasing. This reveals, from the view point of PDF evolution, the center of the wave 
packet may move in the opposite direction of diffusion. As the forwarding of the center the 
width the PDF is always expanding. Therefore, a steady barrier escaping probability can 
always be expected. But the barrier recrossing is greatly reduced. The occurrence of such 


5 






































a nontrivial result may probably be caused by the “Joseph effect” (long-range correlation) 
immersed in the process of fBm {q]. 



FIG. 3: Stationary value Pgt plotted as a function of a. 

For more details, we plot in Figj3]the a dependence of Kgt. From which we can see that K^t 
varies non-monotonically as the increasing of a. There lives a minimum of Kgt at a ~ 0.68 in 
the particular case of what is considered here in this paper. This corresponds to a maximum 
barrier recrossing in the process of fBm which can also be understood from considering on 
the reversed movement of diffusion particles mentioned hereinbefore. Seen from FiglSj both 
{x{t)) and cr^(t) vary slowly as the time goes on. This reveals that the particles will have 
relatively more time wandering in the neighbourhood of the saddle point. Thus results in a 
maximum barrier recrossing and a relatively lower net escaping rate. 

IV. SUMMARY AND DISCUSSION 

In conclusion, we have extended the remarkable Kramers’ rate formula to the case of 
fractional Brownian motion and studied the near barrier dynamics of a reactive fBm pro¬ 
cess. Where several nontrivial properties are revealed such as less barrier recrossing, a 
non-monotonically varying stationary transmission coefficient in great dependent to the 
fractional exponent a and so on. The escaping of fBm particles from the metastable poten¬ 
tial well is found to be more smoothly than the usual cases such as those in one-dimensional 
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non-Ohmic damping systems. And in most cases, the escaping is equivalent to the diffusion 
in the two-dimensional non-Ohmic damping system. We give a brief explanation on the in¬ 
herent mechanism of these results from the view point of PDF evolution. It is expected that 
the present study may have provide more useful information for the in-depth understanding 
of anomalous diffusion dynamics and fBms. 

Yet of course, the courtship concerning on the characteristics of fBm dynamics should be 
more than obtaining a reactive rate. Many facts are still left to be revealed. For an example, 
one can define a reactive flux j{t) = dx{xo = Xb,vo,t)/dt near the saddle point. And find 
the time-dependent varying of it for various a just as we have plotted in FigJU From which 
we can see that in each cases j{t) approaches its maximum quickly and then vanishes to 
zero monotonically in the long time without any other variance. This is reminiscent to the 
remarkable Kramers’ rate problem and can also be understood easily from the view point 
of PDF evolution. 



FIG. 4: Time-dependence of j{t) plotted at various a in which the inserted subgraph is the local 
magnification of it. 

That is to say, at the beginning of the fBm diffusion process, the center of the PDF moves 
quickly to the forward. And at the same time, the width of it expands in a high speed as 
well. Thus results in a rapid emergence of each maximum of j(t). After a short period 
of time, the center of the wave packet may have passed the saddle point. But due to the 
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expansion of its width, the long tail of the wave packet still covers on the barrier. Then a 
small flux j{t) which is close to zero can persistently be found in the long time. Along this 
road, any other dynamical details of the diffusion process can be found without difficulty. 
Therefore we believe that there lives a great possibility for the present study to stimulate 
other in-depth investigations concerning on the fBms in the near future. 
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